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1.  Introduction 


Evaluating  the  positioning  of  assets  such  as  Soldiers  and  equipment  within  a  hostile  environment 
introduces  a  complex  set  of  constraints  defining  the  expected  risk  and  probability  of  success  in  the 
context  of  the  overall  mission  objectives.  These  constraints  generally  require  numerical 
evaluation  and  lack  the  simplicity  of  closed-form  definitions  for  real-world  scenarios.  When  the 
objective  functions  depend  on  the  correlated  positions,  the  result  is  a  complex  geometric 
optimization  problem  that  must  be  carried  out  in  high-dimensional  state  space. 

As  an  example,  we  consider  here  a  prototypical  scenario  in  which  hostile  “red  force”  Soldiers  are 
positioned  within  an  urban  environment,  posing  a  ballistic  (or  sniper)  threat  to  the  intended 
positioning  of  “blue  force”  Soldiers  deployed  to  accomplish  a  prescribed  mission.  The  precise 
threat  must  be  determined  using  a  ray  tracing  technique  to  establish  the  existence  and  distance  of 
a  line-of-sight  (LOS)  path  to  a  prospective  blue  force  location,  with  a  risk  probability  defined  by  a 
ballistic  model  as  a  function  of  distance  (zero  when  occluded).  The  mission  objective  is  further 
defined  as  a  requirement  to  surveil  multiple  targets  within  the  urban  environment.  The  quality  of 
the  observation  from  any  location  in  the  environment  can  be  defined  similarly  to  the  ballistic 
threat  in  that  it  requires  a  LOS  path  and  depends  on  distance  (ignoring  weather  effects  for  now). 
The  optimization  of  the  positioning  of  the  blue  force  Soldiers  must  minimize  the  ballistic  threat 
and  maximize  the  observability  of  targets  for  a  given  scenario.  We  have  developed  a  method  for 
performing  multi-objective  correlated  geometric  optimizations  in  the  context  of  quasi-static 
mission  planning  for  tactical  scenarios. 

This  report  presents  a  formulation  of  a  military  scenario  with  red  and  blue  forces  to  examine  a  ray 
tracing  technique  and  Monte  Carlo  method.  The  background  context  and  the  wording  are 
specific  to  a  military  tactical  operation  for  the  purposes  of  connecting  this  problem  to  a  potential 
real-world  situation.  However,  this  geometric  optimization  method  is  not  limited  only  to  the 
presented  scenario,  but  our  approach  has  the  possibility  to  be  applicable  and  expandable  to  a 
wider  range  of  GIS  intelligence.  Lor  instance,  the  capability  of  computing  timely  urban 
geographical  information  can  be  beneficial  and  effective  for  first  response  teams  in  rescue 
operations  involving  exposure  to  threats  and  desired  targets. 

The  main  contributions  of  this  report  are: 
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•  A  methodology  to  perform  geospatial  analysis  of  a  three-dimensional  (3D)  urban 
representation  to  determine  optimal  geographical  location  given  threat  and  target  positions, 
and 

•  A  Monte  Carlo  method  development  in  the  OpenCL  programming  model  for 
vendor-agnostic  architecture  support  and  future  processor  portability. 


The  remainder  of  this  report  is  organized  as  follows.  Section  2  describes  background  materials 
involved  in  the  implementation  of  the  geometric  optimization.  Section  3  lists  related  work. 
Following  this,  the  methodology  of  geometric  optimization  is  discussed  in  section  4.  The  output 
results  for  the  application  are  revealed  in  section  5.  Finally,  section  6  presents  conclusions  and 
future  work. 


2.  Background 


LOS  calculations  on  a  3D  model  of  an  urban  environment  consisting  of  terrain  and  building 
structures  exhibit  a  characteristic  ripe  for  highly  parallel  processing.  The  triangle  data 
representing  the  geometry  of  a  town  was  extracted  from  a  COLLAborative  Design  Activity 
(COLLADA)  XML-based  file  format,  which  is  used  for  interactive  3D  applications  managed  by 
Khronos  Group  (i).  Casting  only  primary  rays,  this  ray  tracing  (or  casting)  variant  algorithm 
provides  a  held  of  view  from  a  point  selected  in  the  urban  area.  At  its  core,  it  is  the  problem  of 
determining  the  hrst  object  intersected  by  a  ray  as  introduced  by  (4).  3D  LOS  calculations  can 
further  be  extended  to  address  ballistic  hit  probability  and  visual  observation  quality  for  the 
analysis  of  threat  and  exposure.  By  examining  the  computational  geometry  problem  involving 
ray-triangle  intersections  in  a  3D  map,  situational  awareness  can  be  developed  to  augment  the 
geometric  understanding  of  a  region. 

The  parallel  programming  framework  of  OpenCL  (2)  facilitated  the  implementation  of  the 
geometric  optimization  problem  of  determining  the  placement  of  blue  force  Soldiers.  The 
parallelism  of  processing  architectures,  which  is  rampant  and  unavoidable  in  modern  processors, 
is  embraced  from  the  beginning  of  the  development.  As  such,  LOS  calculations  and  the 
optimization  algorithm  are  designed  for  parallel  execution  efficiency.  Details  regarding  the  ray 
tracing  performance,  parallel  programming,  and  graphics  processor  unit  acceleration  are  planned 
to  be  presented  in  a  different  paper.  The  focus  of  this  report  is  to  investigate  the  computation 
methodology  of  a  Monte  Carlo  approach  in  solving  a  geospatial  optimization. 
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3.  Related  Work 


Solving  for  optimal  positioning  of  entities  in  an  urban  environment  falls  under  classes  of 
problems  found  in  computational  geometry.  The  LOS  importance  to  the  problem  allows  for 
some  similarity  to  art  gallery  or  fortress  problems.  Our  problem  shares  the  concern  of  optimal 
placement  in  the  facility  location  problem,  but  with  different  optimization  criteria.  A  review  of 
geometric  optimization  techniques  and  problems  is  presented  in  (3).  To  simplify  the  overall 
objective,  we  refer  to  our  problem  as  one  of  “cooperative  guards.” 

Observations  for  running  Monte  Carlo  schemes  and  Markov  chain  Monte  Carlo  algorithms  in 
parallel  are  provided  in  (ii).  Rosenthal  argues  and  advocates  that  Monte  Carlo  algorithms  are 
ideal  to  parallel  computing.  Lee  et  al.  present  a  case  study  of  running  Monte  Carlo  methods  on 
graphics  cards  that  indicate  potential  for  faster  parallel  simulation  (7).  Our  approach,  however,  is 
to  support  a  Monte  Carlo  method  for  geometric  optimization  on  heterogeneous  platforms  that 
harness  massively  parallel  architectures.  In  addition  to  mapping  for  graphics  processing  units  by 
NVIDIA  and  AMD,  we  envision  targeting  ARM  processors,  Intel  x86  cores.  Digital  Signal 
Processor  (DSP)  architectures,  and  Field  Programmable  Gate  Arrays  (FPGAs). 


4.  Methodology 


To  reiterate,  the  goal  of  the  method  is  to  determine  the  optimum  locations  of  one  or  more 
deployed  blue  force  Soldiers  by  simultaneously  maximizing  the  mission  objective  of  surveilling 
one  or  more  targets  while  minimizing  the  overall  risk  posed  by  one  or  more  enemy  red  force 
Soldiers. 

Let  n  be  the  number  of  blue  Soldiers  used  to  accomplish  the  overall  mission  objective.  Let  the 
position  of  each  blue  Soldier  be  represented  by  a  continuous  position  in  three-dimensional  space, 

rflue)  ^  ^(3)_ 

Let  nired  be  the  number  of  enemy  red  Soldiers  presenting  a  ballistic  threat  to  the  blue  Soldiers. 
Let  the  position  of  each  red  Soldier  be  represented  by  a  continuous  position  in  three-dimensional 
space,  G 
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Let  mtarget  bc  the  number  of  targets  that  are  to  be  observed  as  part  of  the  mission  objeetive.  Let 
the  position  of  each  target  be  represented  by  a  continuous  position  in  three-dimensional  space, 

^{target)  ^  ^(3)^ 

The  objective  function  to  be  optimized  contains  two  components  representing  the  overall  ballistic 
threat  posed  by  the  red  Soldiers  and  the  quality  of  observation  over  all  targets,  both  of  which  are 
functions  of  the  positions  of  the  blue  Soldiers.  The  formulation  here  will  assume  that  the 
locations  of  red  Soldiers  and  targets  are  fixed  during  the  optimization  of  the  blue  Soldier  positions 
(strategic  vice  tactical). 

The  overall  ballistic  threat  probability  is  given  by, 

P(ro,  ri, . . . ,  Tn-i)  =  max  ( maxp^)  (1) 

where  G  [0, 1]  is  the  ballistic  threat  probability  for  blue  Soldier  i  from  red  Soldier  j.  Here 
p  =  0  implies  no  ballistic  threat  vulnerability  and  p  =  1  implies  the  highest  risk  of  ballistic  threat 
vulnerability.  The  rationale  for  the  formulation  in  equation  1  is  that  mission  risk  would  be  based 
upon  the  worst-case  threat  to  any  single  blue  Soldier.  This  formulation  places  a  premium  upon 
the  safety  of  each  individual  Soldier  that  cannot  be  mitigated  by  low  risk  to  the  group  overall.  It 
should  be  noted  that  other  models  could  be  employed  based  on  different  definitions  of  risk. 

The  overall  quality  of  observation  of  all  targets  is  given  by, 

Q{ro,ri,...,rn-i)  =  min  (maxg^j)  (2) 

where  qji  G  [0, 1]  is  the  quality  of  observation  for  blue  Soldier  i  observing  target  j.  Here  g  =  0 
implies  no  observability  and  g  =  1  implies  perfect  observability  sufficient  to  perform  the 
necessary  observation.  The  rationale  for  the  formulation  of  equation  2,  which  differs  in  an 
important  but  subtle  way  from  that  of  equation  1,  is  that  the  overall  quality  of  observation  should 
be  based  on  the  worst  case  for  observing  a  single  target  after  accounting  for  the  best  observation 
overall  from  any  potential  blue  Soldier  observers.  As  an  example,  if  a  single  target  is  left 
unobserved,  this  cannot  be  mitigated  by  good  observability  for  the  remaining  targets.  Note  that 
Qji  and  Q  can  be  interpreted  as  the  probability  of  making  a  successful  observation. 

The  objective  function  that  must  be  minimized  to  determine  an  optimal  solution  is  given  by, 

F  =  WpP  +  Wq{l  -  Q)  (3) 
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where  Wp  and  Wq  are  weighting  coefficients  that  can  be  used  to  set  the  relative  importance  of 
reducing  the  risk  to  ballistic  threat  compared  with  improving  the  observation  of  targets. 

The  above  formulation  puts  no  spatial  constraints  upon  the  locations  of  entities  in 
three-dimensional  space.  In  practice,  realistic  scenarios  will  constrain  the  motion  of  the  blue 
Soldiers  to  a  quasi-two-dimensional  surface  that  can  be  traversed  by  each  Soldier.  This  surface 
can  be  mapped  to  a  two-dimensional  plane  R^‘^\  with  an  elevation  associated  with  each  point  of 
the  plane.  In  order  to  make  the  problem  numerically  tractable,  can  be  discretized  to  form  a 
uniform  two-dimensional  lattice  .  Thus,  the  location  of  each  blue  Soldier  can  be  represented 
as  a  point  Xi  G  since  there  is  a  one-to-one  mapping  from  R^^^  for  a  given  surface  Q. 

In  practice,  there  will  be  constraints  placed  on  the  accessibility  of  locations  within  the  domain. 
Allowable  movement  over  the  surface  can  be  described  by  links  between  adjacent  points  in  the 
lattice  .  Starting  with  a  uniform  and  fully  linked  lattice,  specific  links  are  removed  using  a 
simple  criterion  for  acceptable  movement  of  a  dismounted  Soldier.  If  the  change  in  elevation 
between  adjacent  points  is  greater  than  some  limit,  the  link  between  those  points  is  removed.  By 
constraining  the  allowable  motion  over  this  lattice,  each  point  will  retain  one  to  four  links  to 
adjacent  points.  Based  on  an  initial  location  representing  an  accessible  location  for  a  Soldier, 
only  the  subset  of  points  in  the  lattice  that  can  be  connected  to  the  initial  location  are  retained  as 
the  domain  of  the  problem.  Let  this  domain  be  given  by  C 

As  an  example,  figure  1  shows  the  elevation  map  for  an  urban  environment  and  figure  2  shows  the 
resulting  domain  D^‘^\  where  the  yellow  region  represents  a  constrained  two-dimensional  linked 
lattice.  The  domain  is  constructed  by  starting  with  a  three-dimensional  representation  of  the 
urban  environment  and  choosing  an  initial  location  at  ground  level.  The  connected  lattice 
corresponds  to  all  accessible  locations  from  this  initial  point.  Locations  specifically  excluded  are 
rooftops  and  other  locations  that  would  be  inaccessible  to  a  dismounted  Soldier  without  extreme 
effort. 

The  optimization  problem  reduces  to  determining  the  set  of  positions  Xi  G  that  minimize  the 
objective  function  F  defined  by  equation  3.  Based  on  the  definitions  in  equation  1  and  equation 
2,  this  objective  function  is  not  separable  over  independent  calculations  for  each  blue  Soldier,  but 
is  instead  directly  dependent  upon  the  combined  state  X  =  {xq,  xi,. . . ,  Xn-i}  G  The 

state  space  that  X  belongs  to  is  a  2n-dimensional  space  formed  by  the  n  outer-products  of  D^‘^\ 

In  order  to  determine  an  optimal  value  for  X  in  a  high  dimensional  space,  a  Markov  Chain  Monte 
Carlo  sampling  technique  is  employed.  Using  N  “walkers”  over  the  sample  space,  with  initial 
positions  sampled  randomly  from  a  Markov  Chain  is  constructed  for  each  walker  using 
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N. 


Figure  1 .  Elevation  map  showing  the  input  locations  of  four  enemy 
shooters  located  on  rooftops  (red)  and  four  targets 
located  at  entryways  to  buildings  (green). 


a  standard  Metropolis  algorithm  (9,  5,  5,  6).  A  trial  state  X'  =  {xg,  . . . ,  is  eonstrueted 
from  X  by  drawing  n  random  numbers  Oi  and  updating  eaeh  of  the  state  eomponents  individually 
using  the  following  rule, 


Xi  - 

->■  north 

0<ai<  0.25 

Xi  - 

south 

0.25  <ai<  0.5 

Xi  - 

east 

0.50  <ai<  0.75 

Xi  - 

-r  west 

0.75  <ai<l 

(4) 


where  the  direetions  refer  to  the  links  to  adjaeent  points.  In  the  event  that  a  given  link  no  longer 
exists  within  the  eonstrained  lattiee,  the  aetion  will  return  the  original  point. 


The  objeetive  funetion  F'  is  then  evaluated  for  this  trial  state,  and  the  trial  state  is  aeeepted  based 
on  the  standard  Metropolis  eriterion.  A  random  number  a  is  drawn  from  a  distribution,  and  if 
a  <  F/F',  then  the  trial  state  is  aeeepted  and  X  =  X' .  This  has  the  effeet  of  aeeepting  any  trial 
state  with  a  lower  objeetive  funetion,  and  with  some  probability  aeeepts  trial  funetions  that 
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Figure  2.  Shown  in  yellow  is  the  ground  terrain  that  represents  all 
accessible  points  within  the  scenario  accessible  by 
movement  on  the  ground.  This  region  was  discretized 
and  used  as  a  linked  lattice  for  Monte  Carlo 
minimization. 


actually  increase  the  objective  function  to  avoid  local  minima.  The  optimum  state  is  taken  to  be 
that  whieh  minimizes  F  over  all  states  in  the  Markov  Chain  of  the  N  walkers. 

In  order  to  evaluate  the  objeetive  funetion  for  fixed  red  Soldiers  and  targets,  the  probability  of 
ballistie  threat  and  the  quality  of  observation  is  pre-calculated  for  eaeh  point  within  the 
eonstrained  linked  lattiee  This  proves  to  be  efficient  sinee  it  is  presumed  that  eaeh  point 
within  this  domain  will  be  sampled  at  least  onee  and  each  calculation  is  costly. 

The  ballistie  threat  posed  by  eaeh  of  the  red  Soldiers  is  ealeulated  as  a  probability  field  projeeted 
upon  D^‘^\  There  will  be  rrired  sueh  fields  used  in  the  calculation  of  P.  Equation  1  can  be 
re-expressed  as, 

P(X)  =  max  (  maxp.Txj) )  (5) 

n  \  TUred  / 

where  pj  {xi)  is  the  ballistie  threat  probability  from  red  Soldier  j  at  lattice  point  Xi  G  .  The 
probability  fields  Pj{xi)  are  pre-calculated  using  a  three-dimensional  ray-traced  LOS  method 
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{10),  with  a  polynomial  representation  for  hit  probability  as  a  funetion  of  LOS  distance. 
Examples  of  individual  ballistic  threat  probability  fields  are  shown  in  figure  3. 


Figure  3.  Ballistic  threat  from  each  of  the  four  shooters  as 

determined  by  a  full  three-dimensional  line-of-sight 
ray-tracing  calculation  using  a  mathematical  model  for 
the  ballistic  hit  probability. 


Likewise,  the  quality  of  visual  observation  for  each  target  is  calculated  as  a  probability  field 
projected  upon  There  will  be  mtarget  such  fields  used  in  the  calculation  of  Q.  Equation  2 
can  be  re-expressed  as, 

QiX)  =  min  (maxmfa;,))  (6) 

mtarget  V  n  ^  J'  J 

each  of  the  observation  probability  fields  qi{xj)  is  pre-calculated  using  a  three-dimensional 
ray-traced  LOS,  method  where  the  quality  of  visual  observation  is  represented  as  a  piecewise 
linear  function  of  LOS  distance  with  both  a  minimum  and  maximum  cut-off.  The  minimum 
cut-off  is  introduced  to  prevent  placing  blue  Soldiers  too  close  to  the  targets  since  this  would  not 
make  sense  from  an  operational  perspective.  Examples  of  individual  observability  fields  are 
shown  in  figure  4. 
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Figure  4.  Quality  of  visual  observation  for  each  of  the  targets 

based  on  a  full  three-dimensional  ray-tracing  calculation 
using  a  model  for  the  quality  of  sight  as  a  function  of 
distance. 


5.  Application  and  Results 


The  method  has  been  implemented  using  C-i-i-  and  OpenCL  in  a  code  suitable  for  execution  on 
modern  multi-core  processors  and  many-core  general-purpose  Graphics  Processing  Units 
(GPUs).  Multiple  GPUs  can  be  used  to  evaluate  solutions  at  interactive  speeds,  thus  enabling 
more  tactical  response  times  to  dynamic  situations.  The  details  of  the  implementation  will  be 
presented  elsewhere.  As  a  functional  demonstration  of  the  method,  the  scenario  shown  in 
figure  1  was  constructed  to  include  four  enemy  red  force  Soldiers  placed  on  rooftops  and  four 
targets  (green)  representing  entrance  points  to  a  set  of  buildings.  For  each  red  force  Soldier  and 
target,  the  ballistic  threat  and  observability  fields,  respectively,  were  calculated  using  a  ray-traced 
LOS  algorithm,  with  the  full  three-dimensional  city  map  describing  the  urban  environment. 

These  fields  are  shown  in  figures  3  and  4,  respectively.  The  code  was  then  used  to  perform  a 
static  optimization  to  place  four  blue  Soldiers  in  a  position  to  observe  the  targets  while 
minimizing  the  ballistic  threat  posed  by  enemy  red  Soldiers. 
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A  set  of  4,096  walkers  in  8-dimensional  spaee  (2D  eoordinates  for  four  blue  Soldiers)  were  used 
to  sample  the  ground  plane  linked-lattiee,  eonveyed  in  figure  2,  in  order  to  determine  optimal 
positioning.  A  total  of  409,600,000  samples  were  generated  using  the  eonventional  Metropolis 
algorithm.  Figure  5  illustrates  an  early  state  of  walkers  in  the  Monte  Carlo  sampling,  from  whieh 
shows  sampling  distribution  eoverage  of  the  ground  plane  and  beginning  indieations  of  higher 
density  loeations.  The  final  optimum  positioning  solution,  taken  as  the  state  eonfiguration 
optimizing  the  objeetive  funetion,  is  depieted  in  figure  6  as  the  eolleetion  of  four  blue  Soldier 
loeations. 
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Figure  5.  Early  snapshot  from  a  Markov  Chain  Monte  Carlo 

optimization  taken  after  1,000  steps  using  4096  walkers. 
The  8 -dimensional  state  eoordinate  for  the  4  blue 
Soldiers  is  projeeted  onto  the  2-dimensional  lattiee.  At 
this  early  stage  of  the  optimization,  samples  can  already 
be  seen  to  collect  near  the  locations  that  will  eventually 
be  identified  as  optimal  locations  for  the  four  blue 
Soldiers. 


Performance  of  the  Markov  chain  Monte  Carlo  optmization  OpenCL  kernel  was  collected  for  the 
Intel  Xeon  X5675  and  AMD  Radeon  HD  6970  devices.  For  the  reported  CPU  performance,  the 
execution  time  reflects  a  measurement  for  employing  two  physical  Xeon  CPU  chips  equipped  in 
the  test  workstation.  As  for  the  GPU  results,  execution  times  are  provided  for  one,  two,  and  four 
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Figure  6.  The  solution  of  positions  for  four  reconnaissance 

Soldiers  to  observe  the  targets  (blue).  The  locations  of 
the  reconnaissance  Soldiers  were  determined  using  the 
optimization  based  on  minimizing  ballistic  threat  while 
maximizing  line-of-sight  observation. 


graphics  cards  inside  the  testbed  workstation.  Figure  7  presents  our  performanee  findings  for  the 
geometrie  optimization  algorithm.  The  CPUs  with  12  eores  and  24  threads  at  3.07  GHz  exhibit 
higher  performance  when  eompared  to  a  single  Radeon  card.  Depending  on  an  algorithm’s 
speeifie  eharaeteristies,  optimal  arehiteeture  ean  vary  with  an  elusive  answer.  However, 
maximum  performance  will  depend  on  effeetively  exploiting  the  hardware  parallelism  of 
proeessors.  The  sealing  issue  that  is  observed  with  GPU  eards  ean  be  attributed  to  the  laek  of 
parallelism  for  our  test  ease  of  4,096  walkers.  Not  enough  eomputational  work  is  applied  to  fully 
engage  the  massively  parallel  hardware  design. 
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Markov  chain  Monte  Carlo  Optimization  Performance 


Figure  7.  Execution  time  measured  for  two  Intel  processors  (Xeon 
X5675  six-core)  and  one,  two,  and  four  AMD  graphics 
cards  (Radeon  HD  6970). 


6.  Conclusions  and  Future  Work 


This  report  describes  the  theory  and  demonstrates  the  proof  of  concept  for  calculating  the 
geospatial  optimization  of  objective-constrained  placement  via  Monte  Carlo  method  to  enhance 
tactical  intelligence  in  an  urban  setting.  We  have  presented  a  novel  method  coupling  an 
advanced,  portable,  heterogeneous  parallel  computing  framework  and  a  highly  parallel  ray  tracing 
algorithm,  all  under  the  auspices  of  an  optimization  methodology,  itself,  facilitated  by  massive 
parallelism.  The  end  product  is  an  approach  allowing  for  greater  situational  awareness  through 
advanced  geographical  and  geospatial  processing.  Threat  and  visibility  analysis  in  a  3D  urban 
environment  was  performed  using  a  first-hit  ray  tracing  algorithm.  Once  fields-of-view  are 
generated  for  the  input  entities,  the  numerical  optimization  phase  calculates  the  best  coordinates 
in  the  map  that  minimize  threat  while  maximizing  visibility. 

A  variant  of  a  ray  tracing  algorithm  occupies  the  core  of  the  algorithm.  Here,  the  line-plane 
intersection  computation  was  extended  to  calculate  ballistic  threat  and  quality  of  observation,  but 
this  ray  tracing  approach  can  be  applicable  to  analyzing  radiowave  propagation  in  urban 
operations.  Although  a  specific  scenario  was  presented  in  this  study,  the  core  algorithms  can 
serve  as  a  foundation  to  compute  and  extract  myriad  of  time-critical  intelligence  in  a  3D 
environment. 

To  allow  these  approaches  to  be  used  as  more  than  just  a  before-action  planning  tool  and  work 
more  in  a  real-time  operational  environment,  future  work  also  entails  performing  analysis  on 
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execution  speed  and  examining  differing  processor  architectures  to  achieve  near  real-time 
operation.  With  this  in  mind,  the  source  code  implementation  was  developed  in  OpenCL  to 
conduct  performance  research  on  diverse  parallel  processors  and  accelerators.  Given  the  unclear 
roadmap  for  what  might  become  a  dominant  processing  architecture,  indicated  by  the  rise  of 
ARM  systems  and  accelerators  in  HPC  centers,  our  development  path  decided  to  embrace  parallel 
computing,  vendor  agnostic  portability,  and  heterogeneous  platforms. 

Additional  code  execution  speed  improvements  will  pave  the  way  for  a  more  robust  application. 
Planning  of  extensions  to  dynamic  scenario  optimization  is  underway  where  scenario  assets  are 
not  completely  stationary  within  the  map.  Optimization  will  then  consider  determining  routes 
through  the  map  that  will  satisfy  specific  mission  objectives,  thus  increasing  realism  of  the 
problem.  Determining  the  minimum  cardinality  of  the  set  of  cooperative  guards  is  also  planned 
for  future  exploration.  This  will  be  invaluable  in  determining  not  only  placement  strategies,  but 
also  in  defining  the  minimum  number  of  guards  or  Soldiers  required  based  on  mission 
parameters.  Here  we  will  attempt  to  map  more  completely  the  cooperative  guard  problem  to  the 
fortress  problem,  and  develop  solutions  with  response  times  integrated  with  dynamic  mission 
unfolding  and  streaming  situational  awareness. 


13 


7.  References 


1.  COLLADA  -  3D  Asset  Exchange  Schema,  http  :  /  / www .  khronos  .  org/ collada. 
Accessed:  11/12/2012. 

2.  OpenCL  -  The  Open  Standard  for  Parallel  Programming  of  Heterogeneous  Systems, 
http  :  /  / WWW .  khronos  .  org/ opencl.  Accessed:  11/12/2012. 

3.  Agarwal,  Pankaj  K.;  Sharir,  Micha.  Efficient  Algorithms  for  Geometric  Optimization.  ACM 
Computing  Surveys  (CSUR)  1998,  30  (4),  412-458. 

4.  Appel,  Arthur.  Some  Techniques  for  Shading  Machine  Renderings  of  Solids.  In  Proceedings 
of  the  April  30-May  2,  1968,  spring  joint  computer  conference,  AEIPS  ’68  (Spring),  New 
York,  NY,  pp  37^5,  1968. 

5.  Beichl,  I.;  Sullivan,  E.  The  Metropolis  Algorithm.  Computing  in  Science  Engineering  2000,  2 
(1),  65-69. 

6.  Hastings,  W.  K.  Monte  Carlo  Sampling  Methods  Using  Markov  Chains  and  Their 
Applications.  Biometrika  April  1970,  57  (1),  97-109. 

7.  Lee,  Anthony;  Yau,  Christopher;  Giles,  Michael  B.;  Doucet,  Arnaud;  Holmes,  Christopher  C. 
On  the  Utility  of  Graphics  Cards  to  Perform  Massively  Parallel  Simulation  of  Advanced 
Monte  Carlo  Methods.  Journal  of  Computational  and  Graphical  Statistics  2010,  19(4), 
769-789. 

8.  Liu,  Jun  S.  Monte  Carlo  Strategies  in  Scientific  Computing',  Springer,  corrected  edition, 
January  2008. 

9.  Metropolis,  Nicholas;  Rosenbluth,  Arianna  W.;  Rosenbluth,  Marshall  N.;  Teller,  Augusta  H.; 
Teller,  Edward.  Equation  of  State  Calculations  by  East  Computing  Machines.  Journal  of 
Chemical  Physics  1953,  21,  1087-1092. 

10.  Mdller,  Tomas;  Trumbore,  Ben.  East,  Minimum  Storage  Ray-Triangle  Intersection  J.  Graph. 
Tools  October  1967,  2  (1),  21-28. 

11.  Rosenthal,  Jeffrey  S.  Parallel  Computing  and  Monte  Carlo  Algorithms.  Far  East  Journal  of 
Theoretical  Statistics  2000,  4  (2),  207-236. 


14 


Intentionally  left  blank. 


15 


NO.  OF 

COPIES  ORGANIZATION 

1  DEFENSE  TECHNICAL 
(PDF)  INFORMATION  CTR 

DTIC  OCA 

2  DIRECTOR 

(PDF)  US  ARMY  RESEARCH  LAB 
RDRL  CIO  EL 

IMAL  HRA  MAIL  &  RECORDS  MGMT 

I  GOVT  PRINTG  OFC 
(PDF)  A  MALHOTRA 

1  DIRECTOR 

(PDF)  US  ARMY  RESEARCH  LAB 
RDRL  CIH 
RAJU  NAMBURU 

I  DIRECTOR 

(PDF)  US  ARMY  RESEARCH  LAB 
RDRL  CII 

BARBARA  BROOME 
4  DIRECTOR 

(PDF)  US  ARMY  RESEARCH  LAB 
RDRL  CIH  S 
DAVID  RICHIE 
JAMES  ROSS 
SONG  PARK 
DALE  SHIRES 


16 


